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We analyse the data for the proton structure function F2 over the entire Q 2 domain, including especially low 
Q 2 , in terms of perturbative and non-perturbative QCD contributions. The small distance configurations are 
given by perturbative QCD, while the large distance contributions are given by the vector dominance model and, 
for the higher mass qq states, by the additive quark approach. 
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1. INTRODUCTION 

There now exist high precision deep inelastic 
ep scattering data jl|,|| covering both the low Q 2 
and high Q 2 domains, as well as measurements of 
the photoproduction cross section. The interest- 
ing structure of these measurements, in particular 
the change in the behaviour of the cross section 
with Q 2 at Q 2 ~ 0.2GeV 2 , highlight the impor- 
tance of obtaining a theoretical QCD description 
which smoothly links the non-perturbative and 
perturbative domains. 

In any QCD description of a j*p collision, the 
first step is the conversion of the initial photon 
into a qq pair, which is then followed by the in- 
teraction of the pair with the target proton, see 
Fig.l. Let <t(s,Q 2 ) be the total cross section for 
the process 7*p — * X where Q 2 is the virtuality 
of the photon and yfs is the j*p centre-of-mass 
energy. We can write the dispersion relation in 
the following way: 



cr ( s ' ( 2 2 ) = e, r <5*fe p{s ' m2) w^ 2 ) 



(1) 



where the spectral function p(s, M 2 ) is the den- 
sity of qq states. Following || we may divide 
the integral into two parts: the region M 2 < Qq 
described by the vector meson dominance model 
(VDM) and the region M 2 > Qq described by 
perturbative QCD. To exploit further this idea 
we must achieve a better separation between the 
short and long distance contributions. To do this 
we take a two-dimensional integral over the longi- 



tudinal (z), and transverse momentum (fey) com- 
ponents of the quark, see Figl. 

The contribution coming from the small mass 
region is pure VDM. The part which comes from 
large fcy of the quark can be calculated by per- 
turbative QCD in terms of the known parton dis- 
tributions, whereas for small fcy we will use the 
additive quark model and the impulse approxima- 
tion. That is only one quark interacts with the 
target and the quark-proton cross section is well 
approximated by one third of the proton-proton 
cross section. 

2. The 7*p cross section 

The spectral function p occurring in (Q) may be 
expressed in terms of the 7* — > qq matrix element 
M.. We have p oc |A4| 2 and in terms of the quark 
momentum variables z and k\ the cross section 
in equation ([!]) becomes 

aT =Y, q Q J? Ex = ± 1 / dzd 2 k T (M T M* T )N c ^ImA q - +p 

=E, °i / dzdk T u2+ %C$^ m2q "-w*so ( 2 ) 

3, and e a is 



(<5 +<4> J 

where the number of colours N c 
the charge of the quark in units of e. 

To determine F2(x,Q 2 ) at low Q 2 we have to 
evaluate the contributions to ctt coming from 
the various kinematic domains. First the con- 
tribution from the perturbative domain with 
M 2 > Qq and large fcf., and second from the non- 
perturbative or long-distance domains. 
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Figure 1. The schematic representation of the 
ouble dispersion [l] for the ~f*p total cross section. 
The cut variables M and M' are the invariant 
masses of the incoming and outgoing qq states. 



2.1. The 7*p cross section in the perturba- 
tive domain 

We have to include two graphs, one shown on 
Fig. 2 and the diagonal one, for our calculation 
of the cross section in the perturbative domain. 
Our formula for the cross section has the following 
form: 

ot.l = F T , L (z, k 1T , l T , Q 2 ) ® f{x, l 2 T ) (3) 
where 

f{x,l T )=xdg{x,l 2 T )/d\nl T (4) 

is the unintegrated gluon distribution function. 
The <g> denotes the convolution in the quark 
(z,kir) and gluon ( It ) momenta variables. The 
Ft,l{z,Icit,It,Q 2 ) is the photon-gluon impact 
factor which can be calculated perturbatively. 
From the formal point of view the integrals over 
l T and k T cover the interval to oo. For the l\ 
integration in the domain l T < Iq ~ lGeV 2 we 
may use the approximation 

a s (l 2 T ) f(x,l 2 T ) = f as(l 2 ) f(x,l 2 ). (5) 

The f(x, I 2 ) is the input gluon distribution with 
free parameters which can be adjusted to fit the 
data. 



2.2. Calculating the gluon distribution 

To calculate the perturbative contributions we 
need to know the unintegrated gluon distribution 
f(x, I 2 -). To determine it we carry out the full 
programme described in detail in Ref. Q. We 
solve a "unified" equation for f(x, l T ) which in- 
corporates BFKL and DGLAP evolution on an 
equal footing, and allows the description of both 
small and large x data. To be precise we solve 
a coupled pair of integral equations for the gluon 
and sea quark distributions, as well as allowing 
for the effects of valence quarks. 

Schematically we can write these equations in 
the following form: 

BFKL part DGLAP part 

/ = /° + Ab^T®7+ (Ppg - 1) ® f+Pqg ® S 

S = S°+ S BOX (8) / +P gq ®S (6) 

1<t factorisation 

Following Ref. we appropriately constrain 
the transverse momenta of the emitted gluons 
along the BFKL ladder. There is an indication, 
from comparing the size of the next-to-leading 
ln(l/a:) contribution || to the BFKL intercept 
with the effect due to the kinematic constraint 
H , that the incorporation of the constraint into 
the evolution analysis gives a major part of the 
subleading ln(l/x) corrections. 
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Figure 2. The quark- proton interaction via two 
gluon exchange 
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As in Ref. g we take 1% = 1 GeV 2 , but due to 
the large anomalous dimension of the gluon the 
results are quite insensitive to the choice of Iq in 
the interval 0.8-1.5 GeV. 

The starting distributions for the evolution are 
specified in terms of three parameters TV, A and (3 
of the gluon 



/oOMo) = Nx-\l-x) 



(7) 



where Iq = 1 GeV 2 . The input for the quark sea 
has the following form: 



p _ c-np _ 



2ST P = C. 



(8) 



3. The 7*p cross section in the non- 
perturbative domain 

There are two different non-perturbative con- 
tributions. For M 2 < Q 2 we use the conventional 
vector meson formulae whereas for M 2 > Qq and 
k T < k 2 , we use the additive quark model and the 
impulse approximation. 

3.1. Vector meson dominance part 

As was already mentioned we assume the vec- 
tor meson dominance model to be valid in the 
region where M 2 < Q 2 , and Q\ — 1.5 GeV 2 . We 
include three resonances: p,uj,<j). For complet- 
ness we should also include longitudinal structure 
function Fl{x, Q 2 ). Fl is given by a formula just 
like (0) but with the introduction of an extra fac- 
tor S,Q 2 /My on the right-hand side. £(Q 2 ) is a 
phcnomenological function which should decrease 
with increasing Q 2 . The data for p production in- 
dicate that £( m p) ^;0.7 0, whereas at large Q 2 
the usual properties of deep inelastic scattering 
predict that 



\k\ 



< 



My 



Q 2 



Q 2 



(9) 



Fl 
Ft 

So throughout the whole Q 2 region the contribu- 
tion of Fl is less than that of Ft- In order to 
calculate F L (VDM) we multiply the VDM pre- 
scription for Ft with the factor £Q 2 /M v and use 
an interpolating formula for £ 



My 



My 



Q 



(10) 



with £o = 0.7. 



3.2. Additive quark model and the impulse 
appr oximat ion 

For large produced masses, M 2 > Qq but low 
quark momenta, k T < k$ we use the additive 
quark model and the impulse approximation. Our 
formulae have the following form: 



^ QM =» a / **** ^- +p{ w-) (li) 



where for <Jqq+ p we take, for the light quarks 
2 
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^«+P (W 2 ) = - a p¥ (s = \W 2 ). 



(12) 



To allow for the confinement we replaced Q 
by Q 2 = Q + p 2 in (|ll|), where \i is typically 
the inverse pion radius. We therefore take p 2 = 
0.1 GeV 2 . This change has no effect for Q 2 > p 2 
but for Q 2 < p? it gives some suppression of the 
AQM contribution. 

3.3. The quark mass 

In the perturbative QCD domain we use the 
(small) current quark mass m curr , while for the 
long distance contributions it is more natural to 
use the constituent quark mass Mo- To provide a 
smooth transition between these values (in both 
the AQM and perturbative QCD domains) we 
take the running mass obtained from a QCD- 
motivated model of the spontaneous chiral sym- 
metry breaking in the instanton vacuum S 



A 2 

M ° 2 U 2 +v 



(13) 



The parameter A = 6 1 / 3 /p = 1.09 GeV, where 
p = 1/(0.6 GeV) is the typical size of the in- 
stanton. // is the natural scale of the prob- 
lem, that is p 2 = z(l — z)Q 2 + k\ or p 2 = 
z(l — z)Q 2 + (l T + k T ) 2 as appropriate. For 
constituent and current quark masses we take 
Mo = 0.35 GeV and m curr = for the u and d 
quarks, and Mo = 0.5 GeV and m curr = 0.15 GeV 
for the s quarks. 

4. Summary 

We have made a fit of the cr 7 * p over the entire 
range of Q 2 values. It relies only on the form of 
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the initial gluon distribution see and the bound- 
ary between perturbative and non-perturbative 
contributions. The advantage of our treatment is 
that we use a full set of integro-differential BFKL 
and DGLAP equations. These equations are valid 
over the entire perturbative region. We have 
used very few free parameters which are used to 
parametrise the non-perturbative region. We use 
VDM which is well established in this region. We 
have used running mass prescription in our calcu- 
lation. The growth of m q in the transition region 
is an important non-perturbative effect which we 
find is required by the F 2 data. 

The full and complete calculation of the cross 
section together with the longitudinal part is 
given in |l. 



Fit A: running mass 
Fit A: fixed mass 
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Figure 3. The curves show the virtual photon- 
proton cross section as a function of Q 2 for various 
values of the energy W 
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